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The mechanisms underlying coherent and thermally activated singlet exciton fission in 7r-stacked 
acene crystals are clarified based on quantum dynamics simulations parameterized against a highly 
correlated description of the electronic excitations and their couplings to intramolecular and inter- 
molecular vibrations. In TIPS-pentacene crystals, the relative longitudinal shift of the molecular 
backbones yields large electronic couplings of the triplet exciton pair with both the singlet exciton 
and charge transfer (CT) states. CT-mediated superexchange and direct pathways are found to con¬ 
tribute synergetically to the ultrafast (~100fs) singlet fission process driven by vibronic coherences. 

By contrast, the electronic couplings for singlet fission strictly vanish at the equilibrium 7r-stacking 
of rubrene that exhibits C 2 h symmetry. In this case, the process is incoherent and driven by exci¬ 
tations of symmetry-breaking intermolecular vibrations, rationalizing the experimentally observed 
temperature dependence. 

PACS numbers: 78.47.da, 78.20.Bh, 82.20.Gk 


Singlet excitons in certain molecular crystals such as 
acenes and their derivatives can split into two triplet ex¬ 
citons following a ^in-conserving process known as sin¬ 
glet fission (SF) |ll-[l^. SF has attracted a lot of at¬ 
tention lately in the context of organic photovoltaics, 
since it permits internal quantum efficiency in excess of 
100% through the conversion of high-energy photons into 
two excitons and subsequently two electron-hole pairs at 
donor-acceptor interfaces [j. As triplet excitons usually 
diffuse over longer distances (^/im) compared to singlets 
Si, SF can also be implemented in thick multilayer 
architectures. 

Transient absorption spectroscopy investigations point 
to an ultrafast SF process in molecular crystals of pen- 
tacene i and TIPS-pentacene [l^. Contradicting views 
on the mechanistic aspects of SF in these crystals have 
been reported in the literature Si- Briefly, it has been 
proposed that SF either proceeds through a direct two- 
electron coupling between the singlet exciton (XT) and 
the triplet pair (TT) or follows an indirect mechanism 
where charge-transfer (CT) states act as virtual mediat¬ 
ing states (i.e. superexchange) or are transiently pop¬ 
ulated. In both cases, the singlet-to-triplet conversion 
is believed to be driven by coupling to nuclear degrees 
of freedom and involves either a conical intersection or 
an avoided crossing pathway. Very interestingly, recent 
two-dimensional photon echo experiments performed on 
TIPS-pentacene thin films have demonstrated that: (i) 
while light absorption does not directly generate a co¬ 
herent superposition of XT and TT at time zero, the 
TT population grows in time at the expense of XT with 
an 80fs time constant; and that (ii) yet, the vibrational 
coherences generated upon photoexcitation into XT are 


largely transferred to TT during the SF process [l^ ■ This 
implies an essential role of electron-phonon (vibronic) 
coupling in the coherent SF process. In contrast to TIPS- 
pentacene, SF in the rubrene crystal has been reported 
to be a thermally activated process and is characterized 
by a much longer time constant of few ps @ . The origin 
of the different SF mechanisms in these two families of 
molecular crystals is still unclear. 

In this study, we clarify the SF mechanisms in TIPS- 
pentacene and rubrene based on non-adiabatic quantum 
dynamics simulations fully parameterized against highly 
correlated ah initio electronic structure calculations. We 
find that while ultrafast SF in TIPS-pentacene occurs 
mostly via a CT-mediated mechanism that preserves vi¬ 
bronic coherences, the corresponding process in rubrene 
is incoherent and thermally activated by intermolecular 
symmetry-breaking vibrations (phonons). 

We consider the following linear vibronic coupling 
Hamiltonian in a diabatic representation: 

H = ^/i,(x)|/)(/| + ^ Vij{X){\I){J\ + h.c.){l) 

I i>J 

^/(x) + + (2) 

i i 

where / and J are the indices for the relevant XT, TT, 
and CT states. The diagonal term, hj{x), consists of 
harmonic oscillators for the intra-molecular modes, Xi, 
where uji is the frequency, pi is the momentum, is 
the vibronic coupling, Ej is the excitation energy at the 
ground state geometry. The off-diagonal term, Vij, is 
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FIG. 1: (Color online) (A) Crystal structure of TIPS- 

pentacene and the dimer model for MRMP2 calculations, 
where the side group is simplified as alkyne. The active or¬ 
bitals of dimer are composed of the bonding and anti-bonding 
hybridizations of the HOMO and LUMO of two monomers. 
(B) Crystal structure of rubrene and the dimer model for 
MRMP2 calculations, where the side group is neglected. (C) 
Diagram of electron configurations of XT, TT, and CT states 
in the 3-site, 9-state model. 


the electronic coupling between the states which depends 
on inter-molecular modes, X. The multi-configuration 
time-dependent Hartree (MCTDH) method [l^ has been 
used to solve the Schrbdinger equation, id^ jdt = 77'®', 
with = 'Yhi where T/ are the vibrational wave- 

functions on the respective electronic states. All intra¬ 
molecular vibrational modes and selected inter-molecular 
modes are included in the MCTDH calculations. Our 
model can describe coherent and incoherent SF dynam¬ 
ics mediated by vibronic coupling of realistic systems. 

We consider a molecular trimer, i.e., 3-site, 9-state 
model consisting of 3 XT, 2 TT, and 4 CT states (Fig. 
1C) for the quantum dynamics calculations, where the 
Frenkel exciton on the center site (XT 2 ) is coupled to 
2 TT states similar to the bulk condition. The Hamil¬ 
tonian matrix elements (Table 1) are determined based 
on electronic structure calculations of the dimer mod¬ 
els (Fig. lA and IB) using the multi-reference second 
order perturbation theory (MRMP2) with the correla¬ 
tion consistent polarized basis set (cc-pVDZ), expected 
to provide an accurate description of the excited states. 
The XT, TT and CT energies and their pairwise cou¬ 
plings are obtained by applying a unitary transforma¬ 
tion from the dimer adiabatic electronic states to the 
diabatic representation [T^. The four frontier molec¬ 


ular orbitals, i.e. highest occupied molecular orbital 
(HOMO), HOMO-1, lowest unoccupied molecular orbital 
(LUMO) and LUMO-l-1, are considered as active or¬ 
bitals in MRMP2 calculations on the dimers. The intra¬ 
molecular vibronic couplings are evaluated based on the 
frequencies of normal modes and the geometry optimiza¬ 
tions in the respective states, using the density functional 
theory (DFT) with the PBE functional, where the spec¬ 
tral density is broaden considering the frequency splitting 
[l3| . The GAMESS code [l^ is used for all ab initio cal¬ 
culations. 

First, we address SF in TIPS-pentacene. Table 1 shows 
the diabatic state energies and electronic couplings, as 
obtained from MRMP2 calculations in the ground-state 
geometry. SF in TIPS-pentacene is clearly exothermic. 
Accounting for the intra-molecular reorganization energy 
from the ground state to the exciton {Xxt = 0.15 eV) 
and TT (Att = 0.45 eV) state geometries, the driving 
force AExt-tt amounts to -0.45 eV. 

TABLE I: Hamiltonian matrix (symmetric matrix) of SF in 
diabatic representation where the diagonal and off-diagonal 
elements correspond to the excitation energies and the elec¬ 
tronic couplings (eV), respectively. 
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Ett Ext Ect 
TIPS-Pentacene 1.469 1.618 1.992 
Rubrene 2.451 2.311 3.097 


Vtx Vtc Vxx Vrr Vn Vcc 
TIPS-Pentacene 0.013 0.084 0.018 0.041 0.149 0.017 
Rubrene 0.0 0.0 0.079 -0.175 0.086 0.035 


The slipped-stacked configuration of TIPS-pentacene 
results in large electronic couplings for both the direct 
and CT-mediated pathways (Table 1). Accordingly, the 
quantum dynamics simulations point to an ultrafast SF 
with a timescale of ~100 fs (Fi g. 2 ), in very good agree¬ 
ment with experimental data [l2j| . The contributions 
arising from the two channels can be clearly disentan¬ 
gled by accounting for only one of either the direct two- 
electron XT-TT coupling or the indirect XT-CT and CT- 
TT couplings (Fig. 2B and 2C). Even though the inter¬ 
mediate CT states are ^0.4 eV higher in energy than 
the initial exciton state in our calculations, the superex¬ 
change pathway is found to promote an ultrafast SE pro¬ 
cess and dominates over the direct mechanism, at least 
in the first lOOfs. Such a coherent superexchange can¬ 
not be described by perturbative hopping pictures, high- 
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lighting the importance of including electron-electron and 
electron-vibration couplings on an equal footing when 
solving the time evolution of the system. Note that the 
picture above is not affected by initial excitation condi¬ 
tions, i.e. very similar results are obtained when prepar¬ 
ing the system either in a localized Frenkel exciton or 
in as a delocalized (bright or dark) exciton state (Fig. 
2A, 2D, and 2E). Besides, the ultrafast SF in TIPS- 
pentacene is robust against AExt-tt as long as the 
system is exothermic |14| . 
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FIG. 2: (Color online) Cumulative populations of XT„ 

(blue), CT„ (red), and TT„ (purple) in the quantum dynam¬ 
ics calculations of SF in TIPS-pentacene. (A) Initial exciton is 
localized on the center molecule (XT 2 ), and all the electronic 
couplings are considered. Only one of (B) direct and (C) 
CT-mediated pathways are considered. Dynamics from (D) 
bright delocalized exciton with the initial amplitude, (XTi, 
XT 2 , XT 3 ) = (l/%/3, l/t/3, l/%/3), and from (E) dark exci¬ 
ton, (XTi, XT 2 , XT 3 ) = (1/^3, -1/73, 1/73). 


The coherent nature of SF in TIPS-pentacene is ana¬ 
lyzed based on the overlap of vibronic wave packets on 
the XT and TT states, which provides a measure of the 
coherence M- The substantial vibronic coherence dur¬ 
ing the first ^100 fs (Fig. 3A) indicates that the wave 
packet prepared on the XT potential energy surface is 
only weakly perturbed upon crossing onto the TT hyper¬ 
surface, in line with the coherent nature of the process re¬ 
vealed in the photon echo experiments by Rao et al. [l^ ■ 
Here, the vibronic coupling mediates resonance between 
the initial and final states, and thus purely electronic 
models cannot correctly describe the SF process M- 
The spread of the wave packet on the XT potential sur¬ 
face is wide enough to cross the XT-TT avoided crossing 
immediately after the excitation from the ground state 
(Fig. 3B), hence rationalizing the coherence transfer via 
the direct pathway (Fig. 3A). Our analysis further clari¬ 


fied that the vibronic coherence is transferred throughout 
the superexchange pathway even via the higher-lying CT 
states (Fig. 3A) owing to the strong electronic couplings. 



FIG. 3: (Color online) (A) Overlap integral of vibra¬ 

tional wave packets on the XT and TT states during SF in 
TIPS-pentacene (black), and those considering only the di¬ 
rect (blue) and superexchange (red) pathways, where the ini¬ 
tial exciton is localized on XT 2 . 2Re(tfTT|’FxT) for rubrene 
(green) is negligibly small. (B) Potential energy surface of 
TT (color) and XT (gray scale) states of TIPS-pentacene in 
the effective-modes representation 0 . The solid and dashed 
red circles indicate the Frank-Condon region from the ground 
state and the spread of wave packet, respectively. 

Next, we investigate the origin for the reported ther¬ 
mally activated SF process in rubrene crystals. SF in 
rubrene is exothermic {AExt-tt = - 0.21 eV) when ac¬ 
counting for the reorganization energies {Xxt = 0.18 eV 
and Att = 0.53 eV), although the driving force is smaller 
than pentacene. Remarkably, the TT-XT and TT-CT 
couplings are strictly zero at the C 2/1 equilibrium geom¬ 
etry (Vtx and Vtc in Table 1), hence SF in rubrene 
is expected to be activated by some distortion break¬ 
ing the symmetry. We show below that such an effi¬ 
cient symmetry-breaking mechanism involves relative in- 
termolecular displacements along the rubrene short axis. 
In the adiabatic representation, the C 2 h 7 r-stacking cor¬ 
responds to a conical intersection (zero-gap seam) [T^, 
and the symmetry-breaking induces an avoided-crossing 
between the adiabatic states (Fig. 4C). 

We have thus calculated the potential energy curve 
and the TT-XT and TT-CT electronic couplings with re¬ 
spect to the displacement along the short axis (Fig. 4A 
and 4B), using DFT with the B3LYP functional and the 
Grimme’s dispersion correction. From the linear depen¬ 
dence of the electronic coupling with displacement, one 
can extract the corresponding off-diagonal (non-local) 
electron-phonon coupling, Vjj = XjjX {Xtt-xt = 
0.002 eV, Xtt-ct = 0.003 eV), while the ht with an 
harmonic oscillator of the ground-state potential energy 
curve yields the frequency (~0.007 eV) of the inter- 
molecular symmetry-breaking mode, X. 

The quantum dynamics simulations have been per¬ 
formed by including in the Hamiltonian the symmetry¬ 
breaking inter-molecular mode in addition to the intra¬ 
molecular modes. Temperature effects are built in by 
considering populations of the initial vibrational excita¬ 
tions (i.e., phonon number, Nphonon) of the symmetry¬ 
breaking mode from the Bose-Einstein distribution, 
(Nphonon) = l/{exp(/3tLi) —I}, where /3 is the inverse tern- 
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FIG. 4: (Color online) (A) Potential curve with respect to 
displacement along the short molecular axis as obtained from 
DFT calculations of three rubrene molecules extracted from 
the crystal. To account for the steric hindrance from the other 
side of molecular layer, the mirror image of the calculated po¬ 
tential curve is added. (B) TT-XT (black) and TT-CT (red) 
electonic couplings by MRMP2. (C) Schematic illustration 
of the conical intersection at the C2h 7r-stacking (red circle) 
and the avoided crossing due to symmetry-breaking, where 
the black and blue parabolas illustrate the adiabatic poten¬ 
tials and the seam of diabatic potentials, respectively. (D) TT 
population during SF in rubrene considering a bright initial 
exciton. Two inter-molecular modes (sites 1-2, and sites 2-3) 
are considered, where (Nphonon) are 0 (blue), 3 (black), and 
6 (red). (E) The black line is identical to that in Fig. 4D, 
which is compared with TT evolution considering only one of 
the direct (green) and superexchange (purple) pathways. 


perature. Our quantum dynamics calculations indicate 
an increase in the SF rate with increasing temperature 
(Fig. 4D), consistent with the measured thermally acti¬ 


vated behaviour Q. At 0 K, the vibrational wavefunction 
is localized around the C 2 h region where SF is inefficient, 
since the off-diagonal vibronic couplings {Xtt-xt and 
Xtt-ct) are small in this system. At higher tempera¬ 
tures, the symmetry-breaking mode broadens the wave 
packet over a region where the couplings, Vtx{X) and 
Vtc{X), become large. In contrast to TIPS-pentacene, 
SF in rubrene is: (i) incoherent (Fig. 3A), (ii) mediated 
by the direct 2-electron transition with negligible contri¬ 
bution from the superexchange pathway (Fig. 4E), and 
(iii) occurs around a conical intersection with the wave 
packet on TT being affected by destructive interference 
due to the Berry phase [mil. 

In summary, we have proposed a nonadiabatic quan¬ 
tum dynamical model of SF fully parametrized by first 
principles calculations, and successfully clarified the 
mechanisms underlying the ultrafast and thermally acti¬ 
vated SF in TT-stacked acenes. The present analysis pro¬ 
vides conclusive evidence that in TIPS-pentacene both 
the direct and superexchange pathways mediate SF via 
an avoided-crossing. Notably, the coherent SF via higher- 
lying CT states and its dominant contribution to the ul¬ 
trafast dynamics are not intuitively obvious without the 
explicit quantum dynamical analysis. While the slipped- 
stacked TIPS-pentacene results in strong electronic cou¬ 
plings, the electronic couplings between the TT and other 
states vanish at the equilibrium C 2 h stacking of rubrene. 
SF of rubrene turned out to be driven by thermal exci¬ 
tations of symmetry-breaking vibrations which enhance 
the electronic couplings. The marked difference in the 
SF mechanisms due to the packing symmetry is an in¬ 
teresting feature of 7r-stacked acenes which can inspire 
concepts of molecular design for controlling SF. 
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